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ABSTRACT 

We readdress the calculation of the mass of early-type galaxies using strong gravitational 
lensing and stellar dynamics. Our sample comprises 27 galaxies in the Sloan Lens ACS (SLACS) 
Survey. Comparing the mass estimates from these two independent methods in a Bayesian 
framework, we find evidence of significant line-of-sight mass contamination. Assuming a power- 
law mass distribution, the best fit density profile is given by p oc r - 16 9±o.o5 ^ e s j 10w fa^t 
neglecting the line-of-sight mass contamination produces an overestimate of the mass attributed 
to the lens-galaxy by the lensing method, which introduces a bias in favor of a SIS profile when 
using the joint lensing and dynamic analysis to determine the slope of the density profile. We 
suggest that the line-of-sight contamination could also be important for other astrophysical and 
cosmological uses of joint lensing and dynamical measurements. 

Subject headings: dark matter — galaxies: elliptical and lenticular, cD — galaxies: kinematics and 
dynamics — galaxies: structure — galaxies: fundamental parameters — gravitational lensing 



Introduction 



R ecently the Sloan Lens ACS (SLACS) Sur- 
vey ( Bolton et al.ll2006h obtained images of dozens 



galaxy-scale strong gravitational lenses for previ- 
ously identified early-type galaxies from the Sloan 
Digital Sky Survey (SDSS). The data of both sur- 
veys allow the estimate of the mass of the lens- 
galaxies by two independent methods. 

The mass of a galaxy is one of its most ba- 
sic properties, nevertheless we are only able to 
estimate this quantity based on indirect meth- 
ods. There are several ways to estimate the mass 



methods and their assumptions. 

One of the main motivations of galactic mass 
measurement is to know how it is distributed. 
Some dark matter simulations indicate the emer- 
gence of a universal density profile, which de- 
scribe halos inside a broad mass range of sev- 
eral order s of magnitude, from gal actic to clus- 
ter scales ([Navarro et al.lll996l Il997l) . Simulations 
that include baryons indicate a distinct distribu- 
tion for the stell ar and dark matter components 



of galaxies dFaber & Gallagherl 19791: 


Tvson et al. 




1984; dell' Antonio & Tvson 1996: 


Smith et al. 




2001: Wilson et al. 12001: Guzik & Seliak 2002: 




Prada et al.l2003t Ferreras et al.ll2005: Lintott et al. 


20061 Rettura et al. 20061). and different methods 



insid e a galaxy (jOnorbe et alj l2007t iBaltz et al 
2007). In fact, observations confirm a predomi- 
nance of stellar matter close to the center of galax- 
ies, an d a prevalence of dark matter on the i r out- 



skirts (|Ferreras et al 
Lintott et al.ll2006j) . 



2005- iRettura et all l2006t 



are based on different assumptions and measure- 
ments, therefore it is of great relevance to test 
and compare them. Doing so, we can learn not 
only about the galactic mass, but also about the 



It is of great astrophysical and cosmological in- 
terest to know exactly what is the density profile 
of galaxies, because it can validate or falsify galaxy 
and structure formation models, and even contain 
information about the fundamental pr operties of 
dark matter ( Salvador-Sole et al.ll2007l h 
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In this work we explore the two independent 
galaxy mass estimates obtained using the strong 
gravitational lensing effect and stellar dynamics. 
Both methods are based on very well established 
theories and measurements, and can be combined 
to constrain lens models and the astrophysical 
and cosmological assumptions made. For a the- 
oret ical foundation on t he s ubjects one can con- 
sult Mevlan et al.1 ( 20061 ) and lBinney k. Tremaine 
(|l987h . We are particularly interested in investi- 



gating the role that a possible line-of-sight mass 
contamination can have on the lensing mass es- 
timate, and therefore on the lens-galaxy density 
profile inferred from the joint lensing and dynam- 
ical analysis. 

Work on joint analysis of gravitational lens- 
ing and stellar velocity dispersion in early-type 
l ens galaxies started with few single objects 
Koopmans Treul200ilTreu fc Koopmanskool 



Hamana et al.l 120051 ) . but is rapidly accumulat 
ing events that n ow allow statistical treatment. 
Treu et al. ( 20061 ) studied the distribution of lens 
galaxies in the fund amental plane [also treated by 
Bolton et al.l (|2007j )] and infered that within the 
Einstein radius the isothermal density profile is a 
good approximation for the SL ACS lenses . This 



last re sult was also obtained by IKoopmans et al 



( 20061 ). who quantified the variance in the power- 
law index of the lenses density profile, calculates 
their dark matter fraction and found no signif- 
icant evolution of t he total inner density slope. 
Gavazzi et al.l (|2007l ) used weak gravitational lens- 
ing measurements around lens galaxies that are 
strong lenses to investigate the density profile at 
larger radii, finding again an isothermal mass dis- 
tribution for SLACS lenses. 

Data from SLACS is especially suitable for joint 
strong lensing and dynamic analysis because it al- 
lows precise determination of the Einstein radius 
for each lens-galaxy of a relatively homogeneous 
sample of early- type galaxies. And, at the same 
time, SDSS has precise stellar velocity dispersion 
measurements for all lenses, and redshifts of all 
objects. Previously, the main database for strong 
gravitational lenses containing the same kind of 
data was CASTLE S, the CfA- Arizona Space Tele- 
scope Lens Survey (|Falco et al.lll999l) . however in 
lower number and homogeneity. 

The outline of the paper is the following. In 
Section[2]we present and explain the data used. In 



Section [3] we introduce the theoretical framework 
of mass estimate by strong gravitational lensing 
and by stellar dynamics, and present some models 
for the line-of-sight mass contamination and den- 
sity profile. In Section|3]we lay down the likelihood 
and Bayes analysis formalism used to quantify the 
goodness of fit and compare the various models 
considered. In Section [5] we present our results, 
which are discussed in Section [6] together with our 
conclusions. 

2. The Data 



We use data compiled from IKoopmans et al 



( 20061 ) and Gavazzi et al. ( 2007 ). constructing a 
sample of 27 strong gravitational lensing events 
where the lenses are isolated early-type galaxies 
(E+S0). 

The selected set of galaxi e s is p art of the 
SLAC S Survey (jBolton et all 120061 : iTreu et alJ 
20061 ). which is a Hubble Space Telescope (HST) 
Snapshot imaging survey for strong gravitational 
galactic lenses. The candidates for the HST imag- 
ing were selected spectroscopically from the SDSS 
database and are a sub-sample of the SDSS Lu- 
minous Red Galaxy (LRG) sample. 

For each lens system we are interested in the 
redshift of the background lensed source z s , the 
redshift of the lens z\ , the average stellar velocity 
dispersion inside an aperture a ap and the Einstein 
radius Be- See Table [1] for the compiled data set. 

The source and lens redshifts were determined 
from the SDSS spectra, and the stellar velocity 
dispersion corresponds to the light-weighted aver- 
age inside the 3"diamcter SDSS fiber. 

The Einstein radii were determined from HST 
images using strong lensing modeling of the 
lenses and reconstruction of the unlensed s ources 
(|Koopmans et al.ll2006l : iGavazzi et alJl2007t ). The 
uncertainties on 9e were reported to be around 
5%, so we use this value for all Einstein radii. 
We note that the Einstein radius determined from 
this procedure is a robust attribute of the lens, 
being little sensitive to t he lens model used [see 
Kochanek contribution in iMevlan et all (|2006l )]. 



3. Galaxy Mass Estimate 

We briefly review the independent mass esti- 
mate methods based on strong gravitational lens- 
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ing and stellar dynamics, assuming spherical sym- 
metry. 

3.1. Strong Gravitational Lensing 

When a background source is aligned with an 
intervening mass concentration (lens) the source 
light can be strongly deflected and form a lumi- 
nous ring around the lens (Einstein ring). Less 
well aligned sources can generate arcs and multi- 
ple images. These tracers of the tangential critical 
curve can be used to estimate the mass enclosed 
within the Einstein radius Re = &eDl, 



Mr = 



c z D l D Sq2 



4G D LS 



'E ) 



(1) 



where R>\l,s,ls\ is the angular-diameter distance 
of the lens, source, and between lens and source, 
respectively. 

Note that the source light ray is deflected not 
only by the mass associated to the lens galaxy, 
but by all mass along the line-of-sight, generating 
a contamination that, if ignored, can imply in an 
overestimate of the lens mass. 

We suggest some simple models to subtract the 
light-of-sight mass contamination from the lens 
galaxy mass estimate. We adopt a notation in 
which Ml is the mass estimate from Eq.fTJ) and 
mi, is the lens mass estimate taking into account 
the line-of-sight mass contamination. In the first 
model we subtract a constant contamination from 
the mass estimate of each lens and mj, = Ml — p, 
where p is a free parameter. In the second model 
we assume a constant surface density contamina- 
tion, so, for each lens, the mass overestimate is 
proportional to the lens area, mi — Ml — pA. 
In our third model we propose that the line- 
of-sight mass contamination is not only propor- 
tional to the lens area, but also to the lens mass, 
jxi l — Ml — pMlA. This is an attempt to see 
if more massive lenses have a higher line-of-sight 
mass contamination (this could be the case if they 
are tracing higher density regions). 

The possible dependence of the line-of-sight 
mass contamination with the redshift of the source 
is contemplated by the models m/, = Ml ~ PZl 
and m L = M L — PZlA, where the latter roughly 
mimics the volume covered by the lensed light rays 
from the galaxy lens. 



3.2. Stellar Dynamics 



From Jeans equations (jBinney &: Tremaine 



19871 ) for a spheric halo we have that for any 
isotropic halo profile, the radial stellar velocity 
dispersion is given by 



a 2 (R, z) 



d<$> , 
p-dz, 



(2) 



where we assumed a constant mass to light ratio, 
and <E> is the gravitational potential. 

We define the average line-of-sight velocity dis- 
persion within an aperture of radius R as 



(3) 



where the volumetric mean of the velocity disper- 
sion is calculated inside an infinite cylinder C with 
axis along the line-of-sight, weighted by the den- 
sity profile. 

If we now assume a halo profile of the form 
par 7 and call cr 2 p — (of)(< R ap ) the velocity 
dispersion measured in an aperture of radius R ap , 
then the mass inside the cylinder C [denominator 
in the right side of Eq. fl^)] is 



M D {< R) = -a ap 



2 R 



R 



2+7 



IT- /(7) 



(4) 



where 



/(7) 



1 (5 + 2 7 )(l + 7 )r(-7-l) 

A (3 + 7) r(- 7 -§) 



I)l 2 

2> 



r(-2) 



2 / 

(5) 



Koopmand (|2006t) found similar expressions 

For 7 = —2 (singular isothermal sphere, or 
SIS), Eq.(|4]) reduces to the well known expression 



(< R) 



77 2 D 
Q^apR ■ 



(6) 



Fig. [T] shows the behavior of the factors present in 
Eq. (gl) that depend on 7. For R = R ap the SIS 
profile maximizes the mass contained in a cylinder 
for a fixed average line-of-sight velocity dispersion. 

SLACS data have a velocity dispersion mea- 
sured inside a 3" diameter aperture (9 ap — 1.5") 
and its average Einstein radius is (Oe) — 1.2", 



which imply an average Re/Ru 



0.8. 



Therefore we have two models for the mass en- 
closed inside the Einstein radius estimated from 
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the dynamical method, mo- The first model con- 
siders only a SIS density profile and there is no 
free parameter, tud = M D IS (< Re)- The second 
model allows a free index for the power-law profile, 
m D = M D {< R E ,i)- 

4. Bayesian Analysis 

Since we raise the possibility of various models 
for the dynamic mass estimate and line-of-sight 
contamination, we need a formalism to compare 
these models. So, in addition to a likelihood anal- 
ysis, we also calculate the Bayesian Infor mation 
Crite ria {BIC) and the Bayes evidence (jLiddld 
20071 ). 

The likelihood of both estimate methods to give 
the same galaxy mass is 



U = 



2n(a 2 



: exp 



L,i 



'D.i 



(m L , 



m D , 



2(<7 



2 

LA 



'D,iJ 



(7) 

and the likelihood for a whole set of N galaxies is 
given by 



N 



(8) 



The Bayesian Information Criteria {BIC) is de- 
fined as 



BIC : 



-2 In A 



fclnTV 



(9) 



where C max is the maximum likelihood, k is the 
number of parameters in the model, and N is num- 
ber of data points. 

The Bayesian evidence is defined by 



E = J C(p)P(p)dp 



(10) 



where P(p) is the prior probability distribution 
for parameter p. We use a flat prior probability 
distribution for all the parameters tested. 

To compare two models we calculate the differ- 
ence between their Bayesian Information ABIC 
and A In E. A ABIC > 2 indicates positive evi- 
dence in favor of the model with lower BIC, and 
ABIC > 6 indicates strong evidence. Under ex- 
amination of the Bayesian evidences, AlriE < 1 
does not constitute evidence, 1 < Aln_E < 2.5 
means a significant evidence, 2.5 < Ami? < 5 



means strong evidence, and 5 < Alni5 means de- 
cisi ve evidence in favor of the model with higher 
E (I Jeffreys 111961 . 

The Bayesian evidence is a more rigorou s and 
general criteria than the BIC (jLiddld 120041 ) . but 
we include the latter because it is much simpler to 
calculate, so it is interesting to see if both criteria 
agree. 

We note that Barnabe' fc Koopmansl ( 2007 ) de- 
veloped a framework for joint gravitational lens- 
ing and stellar dynamics analysis using Bayesian 
statistics, but which is very distinct from our ap- 
proach. 

5. Results 

We use the strong gravitational lensing and 
stellar dynamics theories to compute independent 
mass estimates inside the Einstein radius for the 
27 selected galaxies from SLACS. We adopt a con- 
cordance ACDM model (fi m = 0.3 and Q A = 0.7)) 
to provide the redshift-distance relation needed for 
the calculations. 

5.1. SIS Density Profile 

For our simplest model, a SIS density profile 
with no line-of-sight contamination, we find that 
the lensing method overestimates the galaxy mass 
compared to the dynamic method. This is no- 
ticeable in the left plot at FigO A linear fit- 
ting between the two mass estimates gives that 
M L « (1.06 ± 0.08)Mg /s . 

One possible explanation for the observed over- 
estimate of the lensing method in relation to the 
dynamic one is that the former also captures mass 
contaminants in the line-of-sight, while the dy- 
namic method is only sensitive to the galaxy in- 
ternal gravitational potential. We test five mod- 
els for the line-of-sight contamination (models II 
to 15) that have one free parameter and calculate 
their likelihood, BIC and Bayes evidence (see Ta- 
ble [2]). All models with contamination perform 
better than the one with no contamination (model 
10). They have higher likelihood, lower BIC and 
higher Bayesian evidence. 

The line-of-sight mass contamination is signif- 
icant, 4% in model II, 14% in model 12, 11% in 
model 13, 9% in model 14 and 12% in model 15. 
This is computed having as base an average lens- 
ing mass (Ml) = 1.57 x 10 11 M Q , an average lens 
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area (A) = 28. lkpc 2 h 2 and an average lens red- 
shift (zl) = 0.222 for our set of 27 galaxies. 

5.2. Power-Law Density Profile 

If we relax the SIS density profile assumption 
and use instead a power-law density profile for the 
galaxies we obtain models with one more free pa- 
rameter, the power-law index. 

The model with a power-law profile and no line- 
of-sight contamination (model P0) has maximum 
likelihood when its power-law index is —2.02 ± 
0.04, what means that the SIS profile is in fact pre- 
ferred by the system. This occurs because the lens- 
ing mass is larger than the dynamic mass and it is 
an index close to -2, but slightly lower than that, 
which maximizes the dynamic mass inside the Ein- 
stein radius, as is show in FigfT] For our galaxy 
sample {Re /Rap) — 0.8, so the effective curve 
for it is something in between the solid curve for 
R/Rap = 1 and the dotted curve for R/R ap = 0.5. 

All models with power-law profile (models P0 to 
P5) have higher likelihoods than their SIS coun- 
terparts (models 10 to 15, respectively). That is 
expected since the added parameter creates space 
for better fittings, but the extra model complex- 
ity is penalized by the BIC and Bayes evidence 
measures. In fact only for model P2 the gain in 
the likelihood is enough to give a lower BIC and 
higher Bayesian evidence than for its SIC coun- 
terpart, model 12. This model gives the high- 
est likelihood, lowest BIC and highest Bayes evi- 
dence among all models investigated. Therefore it 
is adopted as referential for comparisons with the 
other models. The right plot at FigJ5] illustrates 
how much better is the agreement between the two 
mass estimate methods in model P2. The mass es- 
timate points are closer to the identity line, and 
there is a lower dispersion. A linear fitting to this 
data gives that trl ~ (0.99 ± 0.06)m£>- 

At Figj3] we show the likelihood of model P2 
in its parameter space, the parameter that quan- 
tifies the line-of-sight mass contamination and the 
power-law index of the density profile. There is 
a degeneracy between the two parameters; if we 
take a high power-law index (flatter profile) the 
system demands a higher line-of-sight mass con- 
tamination. In contrast with all other models 
probed, which prefer a SIS profile, model P2 has 
maximum likelihood at a higher power-law index 



7 = —1.69 ± 0.05. As a consequence of this flatter 
profile, farther from the SIS maximization point 
discussed before in relation to FigflJ a higher line- 
of-sight mass contamination is necessary to fill the 
gap between the dynamic mass estimate and the 
lensing estimate. In fact, model P2 has a con- 
siderably higher contaminant surface density than 
model 12 and its line-of-sight contamination rep- 
resents 43% of the average lensing mass. This is a 
very high contamination when compared with the 
ones obtained in the other models, 4% in model 
PI, 14% in models P3 and P4, and 17% in model 
P5. 

6. Discussion and Conclusions 

We used two independent methods, one based 
on strong gravitational lensing and another on 
stellar dynamic, to estimate the mass of 27 early- 
type galaxies from SLACS. Because each method 
has different assumptions and dependencies, the 
comparison of the mass estimate obtained from 
each method enables us to investigate some un- 
derlying aspects of the galaxies and of the lensing 
system, namely the average galaxy profile and the 
line-of-sight mass contamination. 

Likelihood and Bayesian analysis were em- 
ployed to compare models and quantify the evi- 
dence of our results. We find that there is decisive 
evidence of line-of-sight mass contamination for 
the early-type galaxies selected. This is signifi- 
cant, considering that the lens sample was consti- 
tuded by isolated galaxies. Probably because of 
this, we do not find evidence of a higher line-of- 
sight mass contamination for more massive lenses, 
as would be expected if the gal axies belonged to 
group s or c lusters. Previo u sly. iMomcheva et al 



(|2006l ) and IWilliams et~aH (j2006j) also found a 



signi ficant line-of-sight effect on strong lens galax- 



Moustakas et al.l (|2007h found that even in 



underdense local enviroments, the line-of-sight 
contamination may give a considerable contribu- 
tion to galaxy-scale strong lenses, in agreement 
with our results. However, the line-of-sight mass 
contamination is more strongly associated to the 
area of the lens than with its redshift (we also 
tested the dependence with the source redshift 
and lens-source distance and found none). This is 
an indicative that the contamination comes from 
material in the vicinity of the lens and not from 
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along the whole line-of-sight, what is in accordance 
with the cosmological average of zero convergence 
in a random line-of-sight due to the canceling of 
over- and under-dense regions (|Guimaraesll2002t) . 

Our most successful model, which assumes a 
power-law spherical profile and uniform surface 
density line-of-sight contamination, gives evidence 
for a flatter density profile than the one from a 
singular isothermal sphere. In fact, the power-law 
index found is —1.69 ±0.05, but it implies for this 
model a very high line-of-sight mass contamina- 
tion (43%). Such unexpected high contamination 
rises suspition on the limitations of the assump- 
tions and models adopted. 

The main assumptions that we made concern- 
ing the lens galaxies were the sphericity and 
smoothness of the density profiles, no rotational 
support, power-law profile and constant mass-to- 
light ratio along the galaxy radius. These sup- 
positions must be taken as approximations, since 
there is evidence against each one of them, as we 
discuss in the sequence. 

An ellipsoidal symmetry is more general and 
reli able than a spherical one for the len s galax- 
ies (jBaltz et al.ll2007t lHavashi et al.ll2007t ). but it 
is also an approximation and incurs in one ex- 
tra parameter. The existence of substructure on 
galactic scales may also affect the mass estimates 
( Yencho et al.ll2006h . but its inclusio n involves a 
heavy complexltylmore parameters). & et al. 
([20021 ) . assuming spherical symmetry and using 
lensing and dynamic analysis, found a flat inner 
slope for the dark matter distribution in a partic- 
ular galaxy cluste r Thi s result was questioned by 



Meneghetti et al 



( 2005h . who showed the impor- 
tance of ellipticity and substructures using numer- 
ical simulations, 

Contrary to old beliefs, lEmsellem etHI (|2007h 
showed that some early-type galaxies can have 
a significant rotation component, impling some 
rotational support that would require a correc- 
tion of the dynamic mass estimate. The power- 
law approximation for the overall galactic den- 
sity profile has been favored by several works, ap- 
parently as a result of complementary baryonic 
and dark matter profiles, also suggesting a short- 
coming of a constant mass-to-light ratio approxi- 
mation dHamana et al. 2005 ; Ferreras et al 



Lintott et al.ll200a Baltz et al.ll2007f h 



proxi- 
2005; 



Concerning the cosmological framework, the 
model assumed (concordance ACDM) served the 
unique purpose of providing a redshift-distance re- 
lation. If this relation could be obtained indepen- 
dently of cosmic models, then they could be fac- 
tored out of the analysis. Nevertheless, it is the 
very redshift-distance relation that serves as ba- 
sis to existing suggestions of using joint lensing 
and dynamic analysis to probe cosmological mod- 
els and their parameters. In a era of precision 
cosmology, the proper accounting of line-of-sight 
effects would be fundamental in these studies. 

It would be interesting to relax some of our 
assumption and investigate the role of the extra 
parameters. However, models with a large num- 
ber of parameters may provide higher likelihood, 
but lower Bayesian evidence, since the added com- 
plexity must pay its price in a Bayesian sense. A 
larger galaxy sample would probably be necessary 
to reduce the increased degeneracies among the 
free parameters. 

The perspective of more data of the kind nec- 
essary for joint lensing and dynamic studies is en- 
couraging, as well as the potentiality of using them 
to investigate galaxy properties and cosmological 
parameters. On the galactic scale, the density pro- 
file, shape and substructure are the most inter- 
esting targets, and on the cosmological scale the 
redshift-distance relation. Nevertheless, our work 
shows that it is essential for these studies to con- 
sider the line-of-sight mass contamination, since 
neglecting it contribute to an overestimate of the 
galactic mass by the strong lensing method, what 
artificially forces the system to prefer a SIS profile 
because that maximizes the dynamic mass esti- 
mate. 
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Y 



Fig. 1. — The factor present in Eq.Q, which generalizes the mass enclosed inside a radius R for a SIS 
(7 = —2) to a power-law density profile. 
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Fig. 2. — Dynamic and lensing estimates of the mass enclosed by the Einstein radius. The left panel shows 
the results for our simplest model, a SIS with no line-of-sight mass contamination. The right panel is for the 
best power-law density profile model with line-of-sight mass contamination subtraction. 
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Fig. 3. — Likelihood contours in the parameter space of model P2. Likelihood is in units of C max for this 
model, and p is in units of 10 9 M Q kpc~ 2 h 2 . 



Table 1: Compiled data from lKoopmans et al. ( 20061 ) (indicated with a * after the name) and Gavazzi et al 
( 20071 ) (some of which were already present in the former paper) . 



Name 



zi 



0e{") aapjkms- 1 ) 



SDSS 


J002907.8-005550 





227 





931 





82 


228 


± 


18 


SDSS 


J0037-0942 * 





195 





632 


1 


47 


265 


± 


10 


SDSS 


J015758.9-005626 





513 





924 





72 


295 


± 


7 


SDSS 


J021652.5-081345 





332 





523 


1 


15 


333 


± 


23 


SDSS 


J025245. 2+003958 





280 





982 





98 


164 


± 


12 


SDSS 


J033012. 1-002052 





351 


1 


107 


1 


06 


212 


± 


21 


SDSS 


J072805.0+383526 





206 





688 


1 


25 


214 


± 


11 


SDSS 


J0737+3216 * 





322 





581 


1 


03 


310 


± 


15 


SDSS 


J080858.8+470639 





219 


1 


025 


1 


23 


236 


± 


11 


SDSS 


J090315. 2+411609 





430 


1 


065 


1 


13 


223 


± 


27 


SDSS 


J091205.3+002901 





164 





324 


1 


61 


326 


± 


12 


SDSS 


J0956+5100 * 





240 





470 


1 


32 


299 


± 


16 


SDSS 


J095944.1+041017 





126 





535 


1 


00 


197 


± 


13 


SDSS 


J102332.3+423002 





191 





696 


1 


30 


242 


± 


15 


SDSS 


J110308. 2+532228 





158 





735 





84 


196 


± 


12 


SDSS 


J120540.4+491029 





215 





481 


1 


04 


280 


± 


13 


SDSS 


J125028.3+052349 





232 





795 


1 


15 


252 


± 


14 


SDSS 


J1330-0148 * 





081 





711 





85 


178 


± 


9 


SDSS 


J140228. 1+632133 





205 





481 


1 


39 


267 


± 


17 


SDSS 


J142015.9+601915 





063 





535 


1 


04 


205 


± 


43 


SDSS 


J162746.5-005358 





208 





524 


1 


21 


290 


± 


14 


SDSS 


J163028.2+452036 





248 





793 


1 


81 


276 


± 


16 


SDSS 


J223840.2-075456 





137 





713 


1 


20 


198 


± 


11 


SDSS 


J230053. 2+002238 





228 





463 


1 


25 


279 


± 


17 


SDSS 


J230321. 7+142218 





155 





517 


1 


64 


255 


± 


16 


SDSS 


J2321-0939 * 





082 





532 


1 


57 


236 


± 


7 


SDSS 


J234111.6+000019 





186 





807 


1 


28 


207 


± 


13 
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Table 2: Comparison among models and best parameters. 



Model 


m L 




Psup 


Pbest 


Jbest 






ABIC 


Aln£ 


10 


M L 










2.89-10" 


-3U1 


33.7 


14.4 


11 


M l 


-P 


10 11 


(7 ± 2) • 1O 9 M 




2.46-10" 


■299 


28.1 


12.8 


12 


Ml 


-pA 


2-10 9 


(8 ± 2) • 10 8 M e kpc~ 2 h 2 




4.07-10" 


-296 


13.3 


4.1 


CO 


Ml 


— pMlA 


3-10~ 2 


(4.1±0.9)-10~ 3 fcpc~ 2 h 2 




6.95-10" 


-297 


16.9 


6.9 


14 


Ml 


- pz L 


5-10 11 


(7 ± 2) • 10 9 M Q 




5.82-10" 


■298 


21.8 


9.2 


15 


Ml 


- pz L A 


MO 10 


(3.1 ±0.7) • W 9 M Q kpc- 2 h 2 




1.04-10" 


296 


16.0 


5.6 


P0 


Ml 








-2.02 ±0.04 


5.15-10" 


-301 


35.9 


16.2 


PI 


Ml 


-P 


10 11 


(7 ± 3) ■ 1O 9 M 


-1.98 ±0.05 


3.51-10" 


■299 


30.7 


14.4 


P2 


Ml 


-pA 


10 10 


(2.4 ± 0.4) • lO 9 M fcpc" 2 /i 2 


-1.69 ±0.05 


1.61-10" 


-292 








P3 


Ml 


-pM L A 


3-10" 2 


(5± 1) ■ IQ- A kpc~ 2 h 2 


-1.9 ±0.1 


8.97-10" 


297 


19.6 


8.0 


P4 


Ml 


- pz L 


5 ■ 10 11 


(1.0 ±0.3) • IO^Mq 


-1.87 ±0.07 


2.40-10" 


■297 


22.3 


5.1 


P5 


Ml 


— PZlA 


10 10 


(4.2 ±1.0) • 10 9 M e kpc- 2 h 2 


-1.88 ±0.08 


3.46-10" 


296 


16.9 


6.3 



Note. — For models 10-5, mo = M^ IS , and for models PO-5, mo = Mjj(7). The lower limit of the fiat prior for the line-of- 
sight contamination is and the superior limit is p S up, Pbest i s th e best fit value (maximum likelihood). The power-law index 
for the density profile has a flat prior —2.5 < 7 < —1.5 and best fit value -ftest- £max is the maximum likelihood, ABIC the 
difference in the Bayesian Information Criteria, and A In E the Bayesian Evidence ratio, both computed in relation to model 
P2, which has BIC = 1350.3 and E = 5.19 ■ 10~ 295 . 
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